t07Z’3£>b67 



NASA TN 0-6939 



EMISSION OF SOUND 
FROM AXISYMMETRIC TURBULENCE 
CONVECTED BY A MEAN FLOW 
WITH APPLICATION TO JET NOISE 


by Marvin E. Goldstein and Burt M. Rosenbaum 

Lewis Research Center 
Cleveland, Ohio 44135 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION • WASHINGTON, D. C. • SEPTEMBER 1972 



1. Report No. 2. Government Accession No. 

NASA TN D-6939 

3. Recipient's Catalog No. 

4. Title and Subtitle 

EMISSION OF SOUND FROM AXISYMMETRIC TURBULENCE 
CONVECTED BY A MEAN FLOW WITH APPLICATION 
TO JET NOISE 

5. Report Date 

September 1972 

6. Performing Organization Code 

7. Author(s) 

Marvin E. Goldstein and Burt M. Rosenbaum 

8. Performing Organization Report No. 

E-6844 

10. Work Unit No. 
113-31 

9. Performing Organization Name and Address 

Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland* Ohio 44135 

11. Contract or Grant No. 

13. Type of Report and Period Covered 
Technical Note 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, D.C. 20546 

14. Sponsoring Agency Code 


15. Supplementary Notes 


16. Abstract 

A model, based on Lighthill's theory, for predicting aerodynamic noise from a turbulent 
shear flow is developed. This model is a generalization of the one developed by Ribner. 
Unlike Ribner's model, it does not require that the turbulent correlations factor into space - 
and time -dependent parts. It replaces his assumption of isotropic turbulence by the more 
realistic one of axisymmetric turbulence. The implications of the model for jet noise are 
discussed. 



17. Key Words (Suggested by Author(sl) 
Acoustics 
Aerodynamic sound 
Jet noise 

18. Distribution Statement 

Unclassified - unlimited 

19. Security Classif. (of this report) 

Unclassified 

20. Security Classif. (of this page) 

Unclassified 

21. No. of Pages 

41 

22. Price* 

$3.00 


* For sale by the National Technical Information Service, Springfield, Virginia 22151 























CONTENTS 


Page 


SUMMARY . . 1 

INTRODUCTION 1 

BASIC EQUATIONS 2 

TURBULENCE CORRELATION FUNCTIONS 4 

EXPANSION FOR OBSERVATION POINT FAR FROM FLOW 6 

EXPANSION FOR LOW EDDY MACH NUMBER; 

CONVECTIVE AMPLIFICATION 7 

LOCALLY HOMOGENEOUS TURBULENCE; VANISHING 
OF COUPLING TERM 10 

REPRESENTATION OF FOURTH-ORDER CORRELATIONS IN TERMS OF 

SECOND- ORDER CORRELATIONS 13 

AXISYMMETRIC TURBULENCE 15 

AXIS YMMETRIC-QUADRUPOLE CORRELATIONS AND SOUND DIRECTIVITY ... 20 

ISOTROPIC TURBULENCE 26 

AXISYMMETRIC FLOWS 28 

OVERALL INTENSITY 29 

SUMMARY OF APPROXIMATIONS 30 

DISCUSSION OF RESULTS 32 

IMPLICATIONS FOR JET NOISE 33 

CONCLUDING REMARKS 34 

APPENDIX - SYMBOLS 35 

REFERENCES 38 


ni 



EMISSION OF SOUND FROM AXI SYMMETRIC TURBULENCE CONVECTED 


BY A MEAN FLOW WITH APPLICATION TO JET NOISE 
by Marvin E. Goldstein and BurtM. Rosenbaum 
Lewis Research Center 

SUMMARY 

A model, based on Lighthill's theory, for predicting aerodynamic noise from a tur- 
bulent shear flow is developed. This model is a generalization of the one developed by 
Ribner. Unlike Ribner’s model, it does not require that the turbulent correlations factor 
into space- and time -dependent parts. It replaces his assumption of isotropic turbulence 
by the more realistic one of axisymmetric turbulence. The implications of the model for 
jet noise are discussed. 


INTRODUCTION 

In order to use Lighthill's equation (ref. 1) to predict the sound intensity from the 
known properties of a turbulent shear flow, it is necessary to measure or deduce analyt- 
ically the two -point space -time correlations of the second time derivatives of the Reyn- 
olds stresses and then integrate the results over the turbulent region. An evaluation of 
this kind either by exact analytical methods or experimentally in three dimensions would 
be extremely difficult if not impossible. In fact, the only case where the sound emission 
from turbulent flow has been explicitly calculated is for decaying isotropic turbulence 
(ref. 2) which is not directly applicable to jet noise. It is therefore necessary to develop 
approximate models of the turbulence which will allow the noise to be predicted in terms 
of a fairly small number of (possibly experimentally determined) quantities. There have 
been models of this type developed by both Lilley (ref. 3) and Ribner (ref. 4). However, 
both of these treatments are based on the assumption of locally homogeneous and isotropic 
turbulence and Ribner's treatment on the assumption that the space- time correlations of 
the Reynolds stresses factor into the product of a function of space and a function of time. 

In the present report we shall carry through the analysis without making the latter 
assumption. The experimental results of reference 5 show that the assumption of 



isotropic turbulence may not be a good approximation for the correlations of the Reynolds 
stresses needed to predict jet noise. A much better assumption would be that the turbu- 
lence is axisymmetric about the direction of flow. The theory of axisymmetric turbu- 
lence was developed by Batchelor and Chandrasekhar (refs. 6 and 7). In axisymmetric 
turbulence, the statistical properties of the turbulence are the same in every direction 
transverse to the flow. But the properties in the direction of flow can be different from 
those properties in the directions transverse to the flow. The model obtained is more 
complicated than Ribner’s but reduces to his in the case of isotropic turbulence. This re- 
duction occurs without making the assumption that the correlation factor into a space - and 
time -dependent part. This shows that Ribner’s model is actually more general than his 
assumptions indicate. It is shown in the course of the analysis that the decomposition of 
the sound intensity into shear and self -noise with no coupling terms occurs under the rel- 
atively weak assumption that the turbulence is locally homogeneous. (Note that isotropic 
turbulence must also be homogeneous (ref. 8, p. 3). ) In fact, a hierarchy of equations 
obtained by making progressively more restrictive assumptions about the turbulence is 
presented in tabular form. Certain implications of the model for jet noise are discussed. 
It is shown that this model contains both the contributions to the shear noise term postu- 
lated by Lighthill and those postulated by Ribner, thus to some degree unifying these two 
approaches. 


BASIC EQUATIONS 

Lighthill (ref. 1) has shown that at a point x far enough from the flow to be in the 
radiation field of each turbulent eddy (that is, at a distance large compared with (2n)~^ 
times a typical wavelength) the density fluctuations radiated by a localized turbulent flow 
are given by 
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where c Q and p Q are, respectively, the speed of sound and mean density outside the jet 
and the integration is carried out over the region of turbulence. The source function t- 
is given approximately by 
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where v-(y) is the i in component of the total velocity at the point y. Since we are dealing 
with a stationary process in time, we shall suppose (as usual in order to insure conver- 
gence of the Fourier integrals) that r- = 0 for |t| > T where T is some large time 
which will be put equal to infinity at the end of the analysis. We now introduce the 
Fourier transforms 
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Then equation (1) becomes 
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At any point in the far field, the spectral density of the intensity I(x) is 


i 


or 


Pn 2T 


Substituting equation (4) in this expression yields 
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Now upon using the fact that the Fourier transform of a convolution is proportional to the 
product of the Fourier transforms of its components and the fact that the Fourier trans- 
form of f(-t) is the complex conjugate of the Fourier transform of f(t) provided f(t) is 
real, it is easy to see from equations (2) and (3) that 
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where 
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is a fourth -order two -point two-time velocity correlation tensor. 


TURBULENCE CORRELATION FUNCTIONS 


We suppose that the mean velocity is in the yj -direction and changes only in the y 2 - 
direction. * Then 


v^y.t) = 5 u U(y 2 ) + u^y.t) 

where U is the mean velocity and m is the fluctuating (turbulent) velocity. Upon intro- 
ducing this into equation (7), we find that 
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where we have put 


“"(l) “*■” ~ » 

V - -?]' ; = y - y 


(9) 


^i]°k i(y>v> s u ’ 26 ii 6 ij u k u r + u ’ 6 ik 6 n u i u j + u ’ u ' 6 ii 6 ij 6 ik 5 u + u i u j vi 

5 u i u i' u k u ;' - u i u i v'l + 2 u '( 6 ii u j u k u i' - 5 ij u i u k u ;') 


+ u ' u "' 5 ii 6 ik u i u i' + 8 ij 6 u"i u k +26 ij 5 ik u i u i') 


( 10 ) 


^This is a reasonable assumption for a jet if we suppose that y is a locally Carte- 
sian coordinate system at each point of the mixing region. 
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and the double primes indicate that the quantities are to be evaluated at y 1 ' = y' + rj and 
t + r while the primed quantities are evaluated at y* and t. Notice that y'/rj) is 

independent of r. Hence 
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where 8(ce) is the delta function. Since w5(cu) = 0, we see that, when equation (8) is sub- 
stituted into equation (6) and the result substituted into equation (5), the term ^ (y ' , rj) 
makes no contribution to the integral. Hence, upon making this substitution and changing 

the variables of integration from y 1 and y" to y' and ?f for the integration over 

, _ Mf — fl) (2) 

and from y and y to y and for the integration over 3? , we find that 
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But since the variables y" and in the second integral are variables of integration, 
we can replace them by y' and if, respectively. Also since i,j,k, and l are dummy 
indices, we can interchange i with k and j with l without changing the value of this 
integral. Thus, the second integral in equation (12) becomes, after changing the variable 
t to -T, 
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Because the turbulence is assumed to be a stationary process in time, the correla- 
tions must be invariant under a translation in time. Consequently, we can replace t in 
^^^(y',?, -t) by t + r. It is now easy to see that replacing y” by y’,rj-^ by if and 


r by -r in equation (11) merely interchanges the primed and double -primed quantities. 
Hence we see from equations (10) and (11) that 

Upon inserting this relation into equation (13), we find that the second term on the right- 
hand side of equation (12) is merely the complex conjugate of the first. Therefore, 
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EXPANSION FOR OBSERVATION POINT FAR FROM FLOW 


We first introduce the new variable 
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where i denotes the unit vector in the y^-direction. Then 
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Hence, if in equation (14) we put 
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and then change the variables of integration from y 1 and rj to y and r} in performing 
the integration over we obtain 
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where we have put for brevity 
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and the integral over was obtained simply by changing the name of the dummy var- 

variable of integration from y T to y; hence, 
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Because the distance (?7( over which the correlations and ^ijkz are nonzero 

is certainly much smaller than the region of turbulence, we can always suppose that the 
observation point x is sufficiently far away from the flow so that 
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Therefore, upon expanding the integrands in equation (17) and neglecting terms of order 
jvj/jx - yj, we obtain 
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EXPANSION FOR LOW EDDY /V1ACH NUMBER; CONVECTIVE AMPLIFICATION 

We realize that, if the mean flow were zero, the correlation functions 7 (y, rj, r) 
and wou l d de cay so rapidly with rj that we could (as we shall see subse- 

quently) neglect variations in 
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over the range of r\ where 0tt. 
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are nonzero. However, in the case where 


there is a mean flow (see ref. 9), the correlation at the points along the direction of flow 
can extend over long distances since the eddies are moving with the flow and their decay 
time is quite long. In order to compensate for this we introduce a set of coordinates 
which ’’move with the turbulent eddies. ” Thus, we introduce the new variable 


1 = V - U c (y 2 )t 


(20) 


where U c> which is in the y^direction, is the convection velocity of the eddy (ref. 10) 
and following Ffowcs Williams we define the moving axis correlation function by 
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Upon introducing the change of variable (20) into the integral of equation (19) and noting 
that the Jacobian of the transformation is unity, we get 
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where = l^ c !/ c 0 the convection Mach number and 
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is the cosine of the angle between the direction of mean flow and the direction of 
observation. 

Now let l denote a typical correlation length of the turbulence in the moving frame 
and let M g be Mach number characteristic of the turbulent velocity fluctuations. It is 
shown in reference 11 (see sec. 3, eq. (17)) that 
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Now in naturally occurring turbulent flows (i.e. , without artifically induced turbulence) in 
general and jet flows in particular the turbulent velocity is about one -tenth and certainly 
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not more than two -tenths of the mean flow velocity (ref. 12). Consequently, for subsonic 
flows 

— « 1 
c o 

This shows that the exponential in the innermost integral of equation (22) is approximately 
equal to unity for |i;| < l. But by definition t) » 0 for jlfj > l. Hence, we 

can set this exponential equal to unity to obtain the approximate expression 
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Following Ribner (ref. 4) we let l^(xjy) denote the intensity at the point x due to the 
sound emitted from a unit volume at the point y. Then 

z Jx> = f I u<xl y )d y 

and it follows from equation (23) that 

““pq r j r j r k r ; 

32tt 3 c 3 r 6 

where 




-iu!(l-M c cos 0 )t 




r = x - y 

is the vector joining the source point and the observation point and 


is the magnitude of the vector . 

It follows from equations (16), (18), and (21) that 
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the first two integrals on the right-hand side of equation (25) represent the self -noise and 
the shear noise, respectively (refs. 3 and 13). This terminology was introduced by 
Lilley in reference 3 to indicate that the former term represents noise generated by 
turbulent -turbulent interactions whereas the latter term represents noise generated by 
turbulent -shear interactions. The last integral represents a coupling between the shear 
noise and the self -noise. We shall now show, under the relatively mild restriction that 
the turbulence is locally homogeneous and incompressible, that the coupling term in equa- 
tion (25) vanishes. 


LOCALLY HOMOGENEOUS TURBULENCE; VANISHING OF COUPLING TERM 

The theory of incompressible locally homogeneous turbulence is developed in refer- 
ence 8 and in chapter 4 of reference 14. The assumption of local homogeneity implies 
that any two-time correlation function involving the two points y and z, say Q(y,?, r), 
is a function only of rj = z - y and t. Thus 

Q(y,z, t) = q(z - y, r) 


and this implies that 
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On the other hand, 



Hence it follows from equation (20) that 



and by using equation (29), we find that 



We shall now assume that the flow is incompressible and locally homogeneous. 
Then, since 3U/3y^ = 0, the continuity equation is Sm/Sy^ = 0. Hence, if f' is any 
function of the turbulent velocities at z and t + r, 
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which becomes upon taking the time average 



But since the turbulence is locally homogeneous, it follows from equation (30) that 
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Hence 
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Upon using the divergence theorem, this becomes 
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where S is the surface of V and n. is the i in component of the outward unit normal. 
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It follows that, if mf' goes to zero faster than |£|" as |^| — °o, the surface integral 
will vanish in this limit and we obtain 



the integral being carried out over all space. It is shown in reference 15 that 
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Since, by equation (16), rj = i, at r= 0, we can write from equation (28) 

Ri,j k (y,T, t) = ^(iFf 4 ) as |?| — as r=0 (32) 

And since the two observation points always correspond to the same points on the moving 
eddy for all t when is fixed, it follows that the velocities at these points will always 
be at least as well correlated at t = 0 as for any other time delay. Hence we can write 

and it follows from equation (32) that 
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Now we can take f' = u.(y + 77 , t + rju^y + t], t + r) in equation (31) to show that 
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But this shows that the last integral in equation (25) is zero. After changing the names 
of the dummy indices, we can write equation (24) as 
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The first term in this equation is the self -noise term and the second is the shear noise 
term. Ribner (refs. 4 and 13) showed that this decomposition occurs for f I w (x|y)dcu, 
the overall intensity at x, (without including the effect of the moving sources) provided 
the turbulence correlations factor into a time -dependent part and a space -dependent part 
and that the turbulence is locally homogeneous and isotropic. We have shown here that 
this decomposition actually occurs under the less restrictive assumption that the turbu- 
lence is locally homogeneous. 

Equation (33) cannot be further simplified without introducing additional assumptions 
about the turbulence . 


REPRESENTATION OF FOURTH-ORDER CORRELATIONS 
IN TERMS OF SECOND-ORDER CORRELATIONS 

It is argued by Batchelor (ref. 8) that the part of the joint probability distribution of 
the velocities at a fixed time associated with the energy bearing eddies is approximately 
normal at least in so far as the velocities at two points are concerned. This approxima- 
tion is better for some purposes than others. Thus, (ref. 8, p. 176) it gives reasonably 
accurate predictions about the relation between the second- and fourth-order correlations 
This relation is (ref. 8, eq. (8. 3. 11) found to be (see ref. 16, sec. 2. 1. 7 for derivation) 
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But by extending the reasoning used by Batchelor in section 8. 2, we can argue that, when 
the velocity correlations are separated in time as well as in space, their correlation in 
the moving frame will be subject to even more random influences from the neighboring 
flow than when they occur at the same time. In accordance with the central limit the- 
orem, these influences will tend to make the joint probability distribution more normal. 
Hence we expect equation (34), written in the moving frame, to be even more valid when 
r ^ 0 and, in view of equation (30) and (31), we can now write 


R ijkZ^’^’ ^ ~ R ik R jZ + R iZ R jk 
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Upon inserting this into equation (32) and changing the dummy indices, we obtain 
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Having derived this equation in the £ coordinate system, we can now return to the 
rj coordinates and carry out the integration with respect to this variable. We introduce 
the fixed -axis correlation by 
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and the equation can be written in terms of the fixed -axis correlations as 
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In order to simplify equation (36), we now consider a specialization of the turbulence 
field. 


AXI SYMMETRIC TURBULENCE 

Ribner assumed that the turbulence is locally isotropic. However, the turbulence in 

a jet is not necessarily isotropic. For example, the measurements of Jones (ref. 5, see 

fig. 5) show that the correlation R^„(y, 0,0) of the velocities in the direction of flow 

differs noticeably from the nearly equal correlations R 2222^’ 0 ’ 0 ^ and R 3333^’ 0 ’ 0 ^ in 

the directions transverse to the flow (also see ref. 11, sec. 3). These differences could 

have a significant effect on the sound emitted from the jet. A better assumption would be 

that the turbulence is locally axisymmetric. Axisymmetric turbulence was introduced by 

Batchelor in reference 6. He defined axisymmetric turbulence as being any turbulence 

which satisfies the following criteria. "The mean value of any function of the velocities 

2 

and their derivatives must .... be independent of arbitrary rotations of the axes about 
a given direction (which in all practical realizations will be the direction of mean flow) 
and of reflections in planes through that direction. " He later includes the condition that 
it be independent of reflections in planes perpendicular to that direction. 

It follows from Batchelor's definition that any two -point correlation of the velocities 
and their derivatives (with zero time delay) is an axisymmetric tensor (ref. 7). Thus, in 
particular 
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The function should have the same form. 
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Now for any stationary process 


~ZJ~ 


3r n 


u j(y ,t)Uj(y ,t+r) 


B n u. 


-*t=0 


= u.(ir,t) — hy",t) 

-.n 


at 


But by expanding in a Taylor series 
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so that 
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(37) 


is a linear combination of axisymmetric tensors and is, therefore, itself an axisym- 
metric tensor. Thus, when the axis of symmetry is taken in the flow direction (i. e. i the 
T 7 j -direction), the correlation Ry is of the form (ref. 6) 


r) = A 77 . 77 . + Bfijj + C 6 li 6 1 j + D 6 lt ^ + ESj.Tjj 

where A,B,C,D, and E are arbitrary functions of y , r, jrfj = }jvi\ + V 2 + 73 » anc * 
Notice that substituting equation (20) into this equation and recalling that U c is in the 
yj -direction show that the moving axis correlation tensor (27) is also an axisymmetric 
tensor. 

Since the flow is incompressible, 


9Ui(y\t) 


Hence, 


3 


In ref. 6 , Batchelor uses r in place of our 77 and 11 r in place of our tjj. 
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— u i (y’,t)u i (y M ,t + t) = 0 

ayj 1 J 


But since the flow is locally homogeneous, this shows that (see eq. (37)) 

— R ii(y,?f> t) = 0 (38) 

] 

and Rjj is solpnoidal (i.e. , it has zero divergence). 

Since the flow is locally homogeneous, the second -order correlations must be invar- 
iant with respect to translation of the r] coordinate. Consequently, these correlations 

— > A [ f , 1 ^ - T- - , y | » 1 

cannot change if we replace y by y - i 77 : since y = y + - - rj) and y = y + — 

a ^ 2 1 2 

(i, rjj + 77 ), it follows that 

u 1 (y’,t)Uj(y",t + t> - u,(y - it 7, - i - it 1, + i?,t + r) 


Also since the turbulence is stationary in time, the correlations must be independent of 
translations of the time coordinate. Hence, 



ii, I + iy,t + r)-^. 




and in view of equation (36), this shows that 

Rij(y,V, r) = R^(y,-7j,-T) 


We now introduce the axisymmetric tensors 


r) by 



(39) 


(40) 


Then 


where 


R ,i - A m - B \ + c±5 ii 5 ii + *>* v>i - 


n 


2G ± (r) = G(r) ± G(-r) 


for 


G - A,B,C,D, or E (41) 

and the dependence on the other variables besides t has been suppressed for simplicity. 
Then equations (39) and (40) show that 

Ry(y»?, t> = R ji(y, -v, t) (42) 

R[j(y,r],T) = -R”(y, -V, T) (43) 

and equations (38) and (40) show that Rt. and Rj7 are solenoidal. 

It follows by the same argument as used in reference 6 (also see ref. 7, p. 570) that, 
since R^ is an axisymmetric solenoidal tensor which satisfies the symmetry condition 
(42), its coefficients A + , B + , and C + are functions of \ri\ and are even functions of 
and its coefficients D + and E + are equal to one another and are odd functions of 

A similar argument can be used to deduce the dependence of the coefficients of Rf. 
on 77j from the skew-symmetry conditions (43). Thus, by taking i*j and i,j * 1 we 
see from equation (43) that A“(-77 j) = -A - ^). By taking i = j and i,j ± 1, it follows 
that B"(-t?i) = -B”(? 7 ^). By taking i = 1 and j * 1, we find that 

D“(-77 1 ) = E"^) (44) 

Finally by taking i = j = 1, we find that C - ( -? 7 ^) = -C "(r^). 

Now it can be shown by the same argument as in reference 12 that the requirement 
that Rj"j be solenoidal implies that 

Ii®_ + jj + - 1 IJ L. . + 3D” + 77 — + E" = 0 
7] dfj. dr\ 77 dji dr] 

where 77 = |tt| and we have put 77 j = 7771. Since B“ and C" are odd functions of 
(and, therefore, of ju), it is easy to see from this that 

OT-,- SD T-, - 

3D 4 - T] E 

dr] 
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is an even function of 77 Hence, if we put D“ = D '(7?j) - D (- 77 ]) and 
E“ = E~(ri - E - (-77 , it follows that 

3D" 

3D~ o+V ^ + E-=0 

07? 

But it follows from equation (44) that E” = -D~. Hence, 

3D' 

2D o + ^ =0 

37) 

which shows that 

D't? 2 = f(p) 

But as indicated in reference 6, regularity requires that we take f(/i) = 0. Therefore, 
D"(t? 1 ) = D"(-t? 1 ) and equation (44) implies that D"(t)j) - E'^). By adding and subtract- 
ing equation (40) and collecting results, we can conclude that 

Rij (y , V, t) = Rjj (y , 7? , r) + R.'. (y , V, r) (45) 

where 

R ij = A+7 ?i^j + B+6 ij + C + 6 li 6 ij + D+ ( 6 ii ? Jj + (46) 

is an even function of 7 (which follows from eqs. (40) and (41)) whose coefficients A + , 

/ cy rt 

B + , and C + are functions of y, r, and 77 ='\Vi + ^2 + ^3 anc * even ^ unc ^ ions of Vi and 
D + is a function of y, r, and 7? and an odd function of 77 j. Also 

R ij = + Bdjj + + fiytjj) (47) 

is an odd function of r whose coefficients A - , B~, and C” are functions of y, t, and 
77 and odd functions of 77 jj and D“ is a function of y, r, and 77 and an even function 
of 77 j. Notice that in obtaining this representation of the correlation tensor it was only 
necessary to use equations (38) and (39). But it follows from these latter two equations 
and equation (20) that the moving frame correlation tensor must satisfy the same two 
equations with 77 replaced by Hence, this implies that R^ can be expressed in the 
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same form as equations (45) to (47) with £ replaced by r?. 

It can be seen from equations (46) and (47) that RjH = R^. 


AXISYMMETRIC-QUADRUPOLE CORRELATIONS AND SOUND DIRECTIVITY 

We shall now use the symmetry condition derived in the preceding section to simplify 
equation (35). It follows from equation (45) that 


/ 


, -ia)(l-M cos $)t / -Ml-M c °s 0)T|- + 

He / e c R ik R U dr = /?e / e 


/ 


[ R tk R ] + z + R i'k R ^ + Wi + R ik R ti ] d ^ 


Since R^R^ + RT^Rt^ is an odd function of r, we get 


f" e - ioJ 0 " M c COS 6) T (r +r" z + R^R+)d t = -i J sin [w(l - M c cos 6 ) t] + R ik R ]^ dT 

y.oo 

Hence, the left-hand side is purely imaginary and we find 


Re 


f 


-iw(l-M cos 6) r 
e c R 


/*oo 

ik R jf dT = / e 

tX-oO 


-iu>(l-M cos b)t ( , , \ 

K R jW* R ^W 


since the integral on the right-hand side is real. In a similar way it follows that 


-Ml-M. cos 6 )t ( V 2 \ ( V 2 \ / 

^ / e U ( y 2- T ) U (^2 + T )% d - i / to 

(/.OO 


-iu)(l-M cos B)t 


”(’*'t)”( ,j *tK 


dT 


Equation (36) thereby becomes 
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U*|y' ■ 


-ia>(l-M c cos e)' 


( R ik R U + R ik R ' 


w 4 n r.r.r. / -iw(l-M„ cos 6)r f ( V?\ ( ho\ , 

* ■ -TsH 1 / e hv- yH 2 + tK d ” T ( 1 

8ir c Q r J J 

We first consider the first term on the right-hand side of equation (48). It follows 
from equations (46) and (47) that 

R fk R jZ = f [ A \% + B±5 ik + C±6 li 6 lk + D±(6 li ? 7k + 6 lkV] 

X Ja ± t].T] z + + c±6 ij 6 i^ + D± (6i3^ + ( 49 ^ 

Notice that all the products of the coefficients A*, B ± , and so forth which occur in these 
integrands are even functions of ^ an d an d either even or odd functions of rj^. (In 
fact they are all even except those terms which consist of products of D* with either A*, 
B*, or C*. ) But since any integral of an odd function of ^ for i = 1,2,3 is equal to 
zero, it can be seen by inspection that only terms which have two pairs of equal indices 
can contribute to the integral. Thus, the integrals must be of the form 

f R fk R j* ‘ fi i = a uc a ij ( to + b fj 6 ik 6 j; + c tj 6 i; 6 jk < n0 sum on *> i’ k - or l) 


Upon interchanging names of dummy indices, we find that 


VjVj /4 R ,* ^ - 4 r ? r k * b ff?>f - c*rfr , 2 = ( a (j + b fj - < 50 > 


.2 2 


But this shows that the coefficient of r . rf in the summation is just the sum of all terms 

4. 4- — «. J 

in J dr 7 which have two indices equal to i and the other two equal to j. Thus, 

employing the symmetry condition that Ri*; = RiH as noted immediately after equation (47), 

X J J 1 2 2 

we find by inspection of equation (49) that the coefficient of r^r- in equation (50) is 


/( R fj ) 2 ^ * 2 


if i ^ j (no sum on i, j) 
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or 


f( R iif d? ~ 


if i = j (no sum on i) 


Equation (50) now becomes 


: i W f R l R fl \J (R^) 2 d? + r\ J (R± 2 ) 2 d v + (h ^) 2 ^ 

+ 2r l r 3 y*^ 2R 13 R 13 + R ll R 33 )d ^ 

+ 2r l r 2 f^SA + R ll R 22^ dr? 


+ 2r 2 r 3 


f ^ 2R 23 R 23 + R 22 R 33^ ( 51 ) 


Upon carrying out the integration with respect to r\ in terms of polar coordinates, 
we find that 

/fe ) 2 ** = y*( R 3 3) 2 ^ 

f ( 2R 13 R 13 + R 11 R 33 ^ = f ( 2R i 2 R 12 + R ll R 22 )d ^ 
w /'(2R23 R 23 + R 22 R 33^ d ^ = y"( R 22) d7? 

Hence equation (51) becomes 


Vj r k r J f R fk R j7*U r l f( R uf * 2r l( r ! + r :D/" (2R 12 R 12 * R n R 22 )tfi > 


♦ (A ♦ Af /( R±v 


R 22j d7 > 
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2 2 2 2 

and since r j + r 2 + r 3 = r , this becomes 

r i r j r k r ; / B&Bf, d ? = r l y" ( R ll) 2 d ? + 2r l( r2 ' r l) y i 2R 12 R 12 * R ll R 22 ,d ? 

* P - r lf ( R 2 2 ) 2 d ? < 52) 

This can be expressed in terms of the angle 9 between the direction of flow and the di- 
rection r/r from the observation point to the source point (see fig. 1). By using the 
relation 


r i = r cos 9 


(53) 



Figure L - Coordinates for direction between observation and source points. 


equation (52) becomes 


r i r j r k r Z 


/ R fk R j \l d7 ? = cos40 f ( R ll) 2 


d7] + 2 sin 2 0 cos 2 0 J { 2Rj 2 Rj 2 + 


+ sin 4 9 f (r* 2 ) 2 drf 


(54) 


This completes the treatment of the first term on the right-hand side of equation (48). 
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Considering now the second term in equation (48), we find from equation (46) that 
r.r t J U'u"Rt ? cty = rj r t J dy + r 2 J U'U"B + J U’U"C + dy + 2r x r. j U , U ,, D + i Jj dy 

( 55 ) 

where U* = U^y 2 - ^ and U" = |y 2 + The first integral on the right-hand side of 
this equation vanishes whenever odd powers of 77 ^ and 773 occur. Thus, only terms with 
j = l will contribute to this integral. The last integral on the right-hand side of equa- 
tion (55) vanishes unless j = 1 since D + is an odd function of 77 Hence equation (55) 
becomes 

fj r z J uVR* dn = r* J U’u’V^ + r\ j U'u"A + n| drf + r\ J uVA+tj 2 dif 


+ (r 2 + r 2 + 1 - 3 ) J U’U"B + dn + r 2 j u’u”C + d 77 + 2r 2 J uWt/j d rj 


which can be written as 


Yl f u ’ u!,R 37 dn = r \f U , U"Rj 1 d^ + r \ J U'U”R+ 2 dif + r 2 J U’U”R+ 3 (56) 


But it can be seen from figure 1 that 


COS 2 y = — 2 


2 2 
r 2 + r 3 


.2 r 


sin y = 


2 2 
r 2 + r 3 


r 2 + r 3 
sin 2 0 , !lll 3 


and equation (56) can be written as 
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r?r.r 


f UUR^ CI 77 = cos 4 0 f U’U”R 


11 


2 2 2 
+ cos 0 sin 6 cos y 


J U U R ^ 2 cl?] + cos 2 0 sin 2 0 sin 2 


y 


/ U ’ U X3 


d?) 


Using equation (45) and recalling that R f. and R7. are even and odd functions of r 
spectively, we find from equations ( 54 ) and (56) that 


re- 


r i r j r k r Z 


f R ik R jZ d7 ? = bos4e f R u d ^ + 2 sin2 0 cos 2 0 f {, 2R 2 2 + R 11 R 2 2 ) d V 


4 

+ sin o 




+ odd functions of r 


and 


-Li-i y'u’u"R. i d?=cos 4 0 y* u’u"r 11 dfl + COS 2 0 sin 2 0^coa 2 y ^ u’u"R 22 dtl + sin 2 y ^ u'u"R 3g d?j^ + odd functions of 


Recalling that odd functions of r can make no contributions to the integral in equa- 
tion (35), we can how write equation (35) as 


y*|y> - 


16» 3 c*(1 - M„ cos <rtV 

U L* 


.ft 


S®® cos 4 0 + 2 S ®2 cos 2 0 sin 2 0 + s| 2 sin 4 0 ) 


+ 2 js® h cos 4 0 + cos 2 0 sin 2 0 cos 2 y + Sg h sin 2 y 


(57) 


where we have put 


s-=-.^ e f 

t/-QO 


“ la) r T se 
e r -j 6 d t = 




(cod <^ r T)r?j e dr 


(58) 
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(59) 


Sf h 




-iw r 
e 



d t = 



(cos a> r T)r? h dr 


r^ e = J [ R ii R jj + " 5 ip R ij] dT ? (no sum. on i or j) (60) 

r? h s J C U'u"R^ d rj (no sum on i) ( 61 ) 


and where 


o> r = o)(l - M c cos o) (62) 

is the frequency in the reference frame moving with the eddies. Equation (57) expresses 
the spectral density of the far -field intensity emitted from a given point y in a convected 
axisymmetric turbulence. 


ISOTROPIC TURBULENCE 

In the special case of isotropic turbulence, the correlation tensor (37) becomes 

R ij (y , V , t) = A?7.77. + B 6.. (63) 

where A and B are functions only of 77 and t (see ref. 10). By introducing spherical 
cooidinv 3s in equation (60), it can be shown that 

rf- e =r®f for i,j,k, l =1,2,3 (64) 

Upon inserting equation (63) into the continuity condition (38), we obtain 


Hence (ref. 10), if we put 


4A + 77 — 
dr] 


+ i^ = 0 

t] dr] 


Y = B +r] 2 A 


(65) 
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this equation shows that 


A = - ill! (66) 

2 i) St] 

Inserting these equations into equation (63) and then inserting the result into equation (61) 
with i = 2, we find that 



Introducing cylindrical coordinate X, cp, ^ with the axial direction along r \ this 
becomes 



Consequently, 
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It is also easy to see that 


Notice that a similar analysis could be carried through in the moving frame if it were 

assumed that the moving axis correlation tensor (and not the fixed axis correlation ten- 
4 r 

sor ) were isotropic. However, this analysis would lead to the same result. Using 
these results together with equation (64) in equations (57), (58), and (59), we obtain 


yxJy) - 


16tt 3 c 


= (l - M c cos e y 


-,se 

5 11 


2S 1 


sh 4 2 2 2" 

(cos o + cos e sin 0 sin y) 


(67) 


AXI SYMMETRIC FLOWS 

Round jets constitute an important type of configuration for which the flow is axisym- 
metric. In this case equations (57) and (67) which are derived in a locally cartesian co- 
ordinate system at each point y should be averaged over the azimuthal angle y (see 
ref. 4). Since the y-dependence in equations (57) and (67) is explicitly shown, this av- 
eraging is easily accomplished and these equations become 
For axisymmetric turbulence: 


Vxfy) = 


167r 3 c 




M c cos 


cos 4 5 e + 2Sj| cos 2 0 sin 2 e + sin 4 g) 


+ 2 (s® h cos 4 0 + - Sf cos 2 0 sin 2 


( 68 ) 


where 


G sh _ c sh , Q sh 
b T “ b 2 + b 3 

4 

They can not both be isotropic. 

5 

This would actually be more appropriate, since the isotropic tensor must be sym- 
metric in t; and the data presented in refs. 10 and 18 show that the moving axis correla- 
tion tensor is much more symmetric in r than the fixed axis correlation tensor. 
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For isotropic turbulence: 




16tt 3 c 


5 (l - M c cos e) 


4 2 


S 11 + 2S l h ( cos4e + ^cos 2 0 sin 2 e) 


(69) 


OVERALL INTENSITY 

The overall intensity I(xjy) at each field point is obtained by integrating the spectral 
density of the intensity over all frequencies. Thus 


y%00 

!(x|y) = J I w ^Jy)du) 


and equations (57) and (67) to (69) become, respectively, 
General Flow 

Axisymmetrie turbulence: 


I(x|y) = 


8ir 2 c^(l - M c cos 6 ) 5 r 2 


{( a ll cos4e + 2a 12 cos 2 0 sin 2 e + a 22 sin 4 e^ 


[ 4 2 2/ 2 2 Ml 

bj cos 6 + cos 6 sin e (b 2 cos y + b 3 sin y) > (70) 


where 


3 4 se 
a.. = — — r • . 

i] 4 y 

dr 


(71) 


T=0 


bi =iirf h 

' 3x 4 ‘ 


(72) 


T=0 
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Isotropic turbulence: 


I(x|y) - 


8tt 2 Cq (l - M c cos e) 5 r 2 


ja n + 2b 1 (cos 4 0 + cos 2 0 sin 2 0 sin 2 y)J (73) 


Axisymmetric Flow 

Axisymmetric turbulence: 


I(x|y) 


o 


8tt^c^ (l - M c cos e) 5 r 2 


4 2 2 4 \ 

cos 0 + 2a 12 cos 0 sin 0 + a 22 sin 0) 


12 


22 


+ 2 bj cos 4 0 + — (b 2 + bg) cos 2 0 sin 2 0 j (74) 

2 


2 „ 2 , 


Isotropic turbulence: 


l(x\y) = 


87t 2 c^ (l - M c cos e) 5 r 2 




4 „ . 1 2 . 


a^ + 2b t (cos 0 + - cos 0 sin 


sin 2 0^ 


(75) 


SUMMARY OF APPROXIMATIONS 

It is convenient to indicate explicitly the hierarchy of equations which is obtained by 
making the various approximations described previously to obtain the simplified equa- 
tions starting from the essentially exact far -field equation (15). Thus, upon using equa- 
tion (20) and collecting results we arrive at table I. 
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TABLE I. - HIERARCHY OF APPROXIMATIONS 



Isotropic turbulence 
























DISCUSSION OF RESULTS 


/4122\42 

Since 2(cos 0 + - c° s 0 sin 6) = cos 0 + cos 0, equation (75) is exactly the result 
obtained by Ribner (ref. 4, eq. (31)) without including the convective amplification effect 
^1 - M c cos 0)"^- However, Ribner did include the convective amplification effect a 
posteriori (see ref. 4, eq. (34)) by using the higher order correction (in eddy convection 
Mach number) obtained by Ffowcs William in reference 9. It is important to notice that 
we obtained this equation without ever making the assumption that the correlation function 
factors into the product of a space -dependent part and a time -dependent part. Thus, the 
results obtained previously reduce to Ribner ’s results in the special case of isotropic 
turbulence when they are averaged over the azimuthal direction and this occurs without 
making the factorization assumption. 

Now consider the case of an axisymmetric flow. Following Ribner we refer to the 
quantity in braces in equation (74) as the basic directivity pattern of the noise. Since this 
is the directivity pattern with convective amplification and refraction effects neglected, 
the first term 


4 2 2 4 

a^ cos 0 + 2a^2 cos 0 sin 0 + agg sin 0 

4 

is the basic directivity pattern of the self -noise. The coefficient a^ of cos 0 consists 
of the self-correlation of the self -noise part of the longitudinal quadrupole The co- 

efficient ajg depends on the self-c orrelatio ns of the two late ral quadrupoles t ^ anc * 

Tjg and the two cross correlations ( T ii T 22^ S elf anc * ^ T ll T 3Vself' F * na ^y> the coeffi- 
cient a .22 of sin^0 depends on the average of the two self -correlations of the two lon- 
gitudinal quadrupoles T22 and Tgg, on the self-correlations of the lateral quadrupole 
T 2 g, and on the cross correlation of the two longitudinal quadrupoles and the cross cor- 
relation of the longitudinal quadrupoles ^22 T 33' T ** e main contribution comes from the 
self-correlations of the longitudinal quadrupoles. Thus, the basic directivity of the self- 
noise term results from the directivity pattern of a large number of quadrupole correla- 

4 2 2 4 

tions each of which has one of the directivities cos 0, 2 cos 0 sin 0 or sin 0. 

The second term in equation (74) 

4 2 2 

2b j cos 0 + (b 2 + bg) cos 0 sin 0 

4 

is the basic directivity pattern of the shear noise. The coefficient b^ of cos 0 consists 

of the self-correlation of the shear noise part of the longitudinal quadrupole t. ,. The 

2 2 

coefficient (bg + bg) of cos 0 sin 0 consists of the self -correlation of the lateral quad- 
rupoles Tj 2 all d T^g. 
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4 

Thus the cos e term in equation (74) is basically due to a longitudinal quadrupole in 

the direction of flow and the sin^e is basically due to longitudinal quadrupoles trans- 

2 2 

verse to the flow. The cos 9 sin 0 is due to lateral quadrupoles and cross correlations 
of longitudinal quadrupoles. 


IMPLICATIONS FOR JET NOISE 

As we have already indicated, the experiments of Jones (ref. 5) show that the cor- 
relation of the longitudinal quadrupoles in the direction of flow can be much larger than 
the correlation of the longitudinal quadrupoles transverse to the flow direction (see 
ref. 5, fig. 5). Thus, we might anticipate that the cos^e term in the self-noise will be 
considerably larger than the sin^e term. This would cause the self -noise in the mixing 
region to be beamed downstream. Since the spectrum of the shear noise peaks about one 
octave below the self -noise, t'his would at first seem to contradict the fact that the basic 
directivity of the high frequency sound tends to be nondirectional whereas the low fre- 
quency sound tends to be beamed downstream. However, this contradiction might be, re- 
solved when one considers the spectrum of the self -noise term. The small scale eddies 
will contribute mostly to the high frequency part of the spectrum of the self -noise. But’ 
these eddies tend to be much more isotropic than the large scale eddies which contribute 
mostly to the low frequency part of the spectrum. Thus the high frequency part of the 
self-noise spectrum will still be isotropic and only the low frequency part will be non- 
isotropic. This nonisotropy will result in a downstream beaming of the low frequency 
part of the self -noise so that both the low frequency part of the self -noise and the shear 
noise (which is all low frequency) will be beamed downstream. However, the high fre- 
quency part of the self -noise spectrum will be nondirectional. Hence, just as in the case 
of isotropic turbulence, we find that the high frequency sound is nondirectional while the 
low frequency sound is beamed downstream even though the self -noise term itself may be 
highly directional. The model of isotropic turbulence implies that this term is basically 
nondirectional (see eq. (75)). 

We expect that the effects of nonisotropy of the turbulence can be significant in de- 
termining the directivity. pattern of the narrow band noise as anticipated by Ribner 
(ref. 4, sec. 9). 

For Ribner 's model with isotropic turbulence (ref. 4) only the self -correlations of 
the longitudinal quadrupole and the lateral quadrupole contribute to the shear 
noise. In contrast, Lighthill (ref. 17) argued that the shear noise should be due almost 
entirely to the self -correlation of the lateral quadrupole which does not occur in 
Ribner 's theory. However, the fuller pattern predicted by Ribner 's result is in closer 
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agreement with the experimental data (ref. 18). We see that, for axisymmetric turbu- 
lence, all three quadrupoles are present. This unifies these divergent viewpoints to 
some degree. 


CONCLUDING REMARKS 

A model has been developed to include the effects of nonisotropic turbulence in pre- 
dicting jet noise. The turbulence model chosen is that of axisymmetric turbulence about 
the direction of flow. Thus, the statistical properties of the turbulence are independent 
of rotations about the direction of flow. The model reduces to Ribner's model in the spe- 
cial case when the turbulence is isotropic. However, Ribner’s result is recovered with- 
out assuming, as Ribner did, that the turbulence correlation functions factor into a space- 
dependent and a time -dependent part. This shows that Ribner's model is actually more 
general than the assumptions he used to derive it would indicate. It is also shown that the 
only assumption required to assure that the intensity will be equal to the sum of a shear 
noise term and a self -noise term (with no cross coupling terms) is the relatively weak 
one that the turbulence be locally homogeneous. 

It is indicated that the basic self -noise pattern may actually be more highly direc- 
tional than previously anticipated and this can occur without contradicting the observed 
result that the low frequency noise is beamed downstream. The shear noise terms for 
axisymmetric turbulence contain both the quadrupoles of Ribner's model and those pro- 
posed by Lighthill which are missing from Ribner's model. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, May 26, 1972, 

113-31. 
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APPENDIX - SYMBOLS 


A 

A ± 

a ij 

B 

B* - 

b i 

C 

c* 

c o 

D 

D ± 

E 

E ± 

F 

G, G* 

I(x) 

I(x|y) 

l 

M c 

M e 

Q 


q 


R ijkZ ,R i, jk’ R ij 



» 

3t 


ijW» 
( 1 ) 
ijk V 


f ijkZ 

a(2) 



coefficient in axisymmetric tensor 
even/odd part of A in t 
fourth time derivative of correlation function 
coefficient in axisymmetric tensor 
even/odd part of B in r 

fourth time derivative of correlation function for shear noise 
coefficient in axisymmetric tensor 
even/odd part of C in t 
speed of sound 

coefficient in axisymmetric tensor 
even/odd part of C in t 
coefficient in axisymmetric tensor 
even/odd part of D in r 
defined in eq. (65) 

any coefficients of axisymmetric tensor (see eq. (41)) 

intensity at point x of sound field 

intensity at point x of sound emitted from point y 

spectral density of I(x) 

spectral density of I(xjy) 

correlation length of eddy 

convection Mach number 

eddy Mach number 

arbitrary function of y , z, and t 

arbitrary function of z-y and r 

moving frame correlation functions of turbulent velocities 
fixed frame correlation tensor of turbulent velocities 
turbulence correlations in moving frame 
turbulence correlations in fixed frame 



r 

r 


r i 

r f 

se 

rr. 

i] 

Sf h 


S. 


se 

ij 


T ij 

t 

U 

U_ 


u 


u i 


X 


x i 


y 


y ,y 

y i 

z 


y 

A 

? 

0 

P 

£*o 


x-y, vector connecting source and observation points 

magnitude of r 

component of r 

shear noise correlations 

self -noise correlations 


sh 


proportional to Fourier transform of r^ 

cp 

proportional to Fourier transform of r.. 
Fourier transform of velocity correlations 
time 

mean velocity 
convection velocity 
turbulent velocity 
component of u 
component of total velocity 
coordinate of observation point 
component of x 


( » y 2 

average source coordinate; y = ly^ _ 
coordinates of soirrce points 

ii. 

i 1 11 component of y 
coordinate of source point 
component of z 
azimuthal angle 

Fourier transform of density fluctuation 
^-U c r 

correlation distance y” - y 1 


y 2 + y y 3 + y 3 


-V 

component of rf 

angle between observation point and direction of mean flow 
instantaneous density 
mean density 
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t time delay 
c j j frequency 


Superscripts: 

denotes quantity evaluated at point y* and time t 

denotes quantity evaluated at point y" = y’ + r\ and time t + r 
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